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Abstract 

Although the standard generally-covariant Dirac equation is unique 
in a topologically simple spacetime, it has been shown that it leads to 
non-uniqueness problems for the Hamiltonian and energy operators, 
including the non-uniqueness of the energy spectrum. These problems 
should be solved by restricting the choice of the Dirac gamma field in 
a consistent way. Recently, we proposed to impose the value of the 
rotation rate of the tetrad field. This is not necessarily easy to im- 
plement and works only in a given reference frame. Here, we propose 
that the gamma field should change only by constant gauge transfor- 
mations. To get that situation, we are naturally led to assume that 
the metric can be put in a space-isotropic diagonal form. When this 
is the case, it distinguishes a preferred reference frame. We show that 
by defining the gamma field from the "diagonal tetrad" in a chart in 
which the metric has that form, the uniqueness problems are solved 
at once for all reference frames. We discuss the physical relevance of 
the metric considered and our restriction to first-quantized theory. 



1 Introduction and summary 

Dirac's original equation, which describes the behaviour of elementary spin 
half particles, is valid only in Cartesian coordinates on a Minkowski space- 
time. It thus has to be modified in order to be written in a curved spacetime, 
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or already in a flat spacetime with non-Cartesian coordinates, e.g. in a ro- 
tating frame, as is the case of our Earth. This leads to redefining the Dirac 
matrices 7 M as a field depending on the position X in the spacetime V. That 
field is a priori only assigned to satisfy the basic anticommutation relation 



Here, (g^) is the inverse matrix of the matrix (g^u) made with the compo- 
nents g^v of the metric g in a local coordinate system or chart \ : X t— > (a; M ). 
With respect to the original (Dirac's) Dirac equation, we have thus the 
new possibility that now different admissible fields of Dirac matrices, say 7^ 
and 7 M , are related together through a point- dependent invertible complex 
matrix, which is unique up to a complex scalar \(X) [1, 2j: 



The point-dependent similarity transformation S [3J is usually called a "lo- 
cal" similarity transformation (e.g. [1]). It should not alter the physical 
properties of the fermions described by the modified Dirac equation, and is 
thus a local gauge transformation. The standard modification of the Dirac 
equation for a curved spacetime is due to Fock and to Weyl; we call it the 
Dirac-Fock- Weyl (DFW) equation so as to distinguish it from both the orig- 
inal Dirac equation and the alternative modifications of the latter, which 
were proposed in previous works [61 [7J. For the DFW equation, the local 
gauge transformation S is deduced from a local Lorentz transformation ap- 
plied to a tetrad field. This implies that, for DFW, S(X) at each spacetime 
point X G V belongs to the spin group Spin(l, 3). The DFW equation is co- 
variant under any such transformation, provided that it be smooth J3j [5j [8] . 
(This is in addition to being covariant under a general coordinate change, or 
"generally-covariant".) Thus, the DFW equation is unique in a given curved 
spacetime, at least [5] if the latter is topologically simple. 



1 The Dirac matrices 7^ depend on the chart as is obvious from Eq. ([IJ , and in general 
a chart is defined only locally, i.e. in some open subset U of the spacetime V. However, the 
7 P matrices are made with the components of the chart-independent "7 field" . The latter 
is a global section of the vector bundle TV <S) E E°, with E the spinor bundle (Sect. [2]) 
and E° its dual vector bundle [6]. Using the coordinate basis (<9 M ) associated with a chart, 
a local frame field (e a ) on E, and its dual frame field (9 a ) on E°, the local expression of 
the 7 field is 7 = 7^° <9 M <£> e a ® 9 b ■ The Dirac matrices are 7^ = (7^°)- Considering the 
7 field leads one immediately to the correct transformation law of the 7^ field [H] . 



7Y + 7Y = 2ri 4 , /x,z/G{0,...,3}. 



(1) 



T» = S- x <fS, /i = 0,...,3. 



(2) 
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Nevertheless, it has been recently discovered [5] that, in any given refer- 
ence frame (see below), or even in any given coordinate system, neither the 
Hamiltonian operator H nor the energy operator E is unique for the DFW 
equation, and instead both H and E depend on the choice of the tetrad field, 
in such a way that even the energy spectrum is not unique. It had been 
previously observed on a particular case by Ryder [TU] that the presence or 
absence of Mashhoon's "spin-rotation coupling" term [TTJ [T2] in the DFW 
Hamiltonian may depend on the choice of the tetrad field. Thus, the gauge 
invariance of the DFW Lagrangian, that is its invariance under the spin 
transformations S : V — > Spin(l,3), is not a sufficient protection against an 
unphysical dependence (e.g., of the energy spectrum) on the arbitrary choice 
between equally valid fields of Dirac matrices. To solve this non-uniqueness 
problem, we have to restrict the choice of the tetrad field. The generic non- 
uniqueness of the Dirac Hamiltonian and energy operators is true also j9] 
for the alternative generally-covariant Dirac equations [6j [7J, for which the 
7 M field need not be deduced from a tetrad field. The standard "Dirac La- 
grangian" valid for DFW can be extended to these alternative equations, see 
Eq. (|74|) below, and is invariant under those local gauge transformations 
S : V — > GL(4, C) which satisfy a certain PDE [9j [IB] . However, in the 
present work, we define the 7^ field from a tetrad field, not only for DFW 
but also for that alternative equation which we consider. Thus, the local 
gauge transformations are restricted to Spin(l,3) — see Eq. (f4"3"|) below. 

In a recent work [13] , a solution of this non-uniqueness problem has been 
proposed, in particular for DFW, based on the following observations. The 
Hamiltonian and energy operators are defined with respect to a given refer- 
ence frame [16J — the latter notion being formally defined as an equivalence 
class F of charts, any two of which exchange by a purely spatial coordinate 
change [17J. The data of a reference frame F fixes a unique four- velocity 

2 The work [9] discusses this problem in detail. The work |13) begins with a summary of 
the problem and its cause, and analyzes an attempt [141 115] at solving the non-uniqueness 
problem of the DFW Hamiltonian operator: in a given coordinate chart, there are different 
"Schwinger tetrads" , related two-by-two by, in general, a time-dependent rotation. {This 
can be seen already for a Cartesian chart in a Minkowski spacetime, see Eq. (96) in 
[T5] — Eq. (94) in the arXiv version.} That gives by Eq. (15) in [13 a time-dependent 
gauge transformation, whence two physically non-equivalent Hamiltonians, Eq. (17) in 
[13j . Moreover, the non- uniqueness generally subsists even if a unique tetrad field is fixed 
in a given chart: see Sect. |4] below, and see point (iii) in Sect. [6l 
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field v [TSJ [T7], and also fixes a unique rotation rate field O [T5J [TH1 [Sl- 
it is natural to impose on the tetrad field (u a ) (a = 0, ...,3) the condition 
that the time-like vector of the tetrad is the four-velocity of the reference 
frame: Uo = v [201 1211 [13] . It was shown [13] that then the spatial triad (u p ) 
(p = 1, 2, 3) can only be rotating w.r.t. the reference frame. All tetrad fields 
(u a ), which have the same time-like vector field uo, and for which the rotation 
rate tensor field 3 of the spatial triad (u p ) is the same, give rise to equivalent 
Hamiltonian operators as well as to equivalent energy operators |13] . Two 
natural ways to fix the tensor field 3 are: i) S = ft, where ft is the unique 
rotation rate field of the given reference frame, and ii) 3 = 0. Either choice, 
i) or ii), thus provides a solution to the non- uniqueness problem. These two 
solutions are not equivalent, so that experiments would be required to de- 
cide between the two. The difference between the two solutions is likely to 
be of the order of the Mashhoon spin-rotation coupling term [T0| [13] , which 
currently is still too small to be tested. It may also be difficult to implement 
these solutions in practice, especially the solution i) which involves the tensor 
field f2. Moreover, each solution is valid only in a given reference frame |13j . 

The aim of the present paper is to present a simpler solution of the non- 
uniqueness problem of the generally- covariant Dirac theory. This solution 
consists in restricting the choice of the tetrad field to a class, inside which 
any two choices are related together by a constant Lorentz transformation of 
the first tetrad field, which translates into a constant gauge transformation 
S of the 7^ field. This is the closest possible solution to the situation for the 
original Dirac equation in a Minkowski spacetime, in which any two sets of 
Dirac matrices exchange by a constant similarity transformation [TJ |2] . The 
idea of restricting oneself to 7^ fields which exchange by constant similarity 
transformations is quite obvious from special relativity, but its application 
to the general Dirac equation in a curved spacetime is not. First, in Sec- 
tion [2] we show from previous work [9] that this idea can work a priori only 
for two versions of the generally-covariant Dirac equation: DFW and the 
so-called "QRD-0" version [6j. For each of these two versions, two 7 M fields 
that exchange by a constant gauge transformation give rise to equivalent 
Hamiltonian operators, as well as to equivalent energy operators, and this in 
any possible reference frame. Then in Section [3] we present a simple general 
prescription to select a tetrad field in a given reference frame. We show in 
Section H] that in general this prescription leaves us with the non- uniqueness 
problem, even if the metric is diagonal — in which case this prescription coin- 
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cides with the widely used "diagonal tetrad" prescription §M§, which is thus 
non-unique in general. To get it unique leads one naturally, as we show, to 
consider a space-is otropic diagonal metric: 

(<7 M „) = diag(/, -h, -h, -h), f>0,h>0. (3) 

For that metric, we prove in Section [5] that the "diagonal tetrad" prescription 
is unique. [f| In fact, that prescription applied to that metric leads exactly to 
the desired property. Namely, if one applies the "diagonal tetrad" prescrip- 
tion (I24p in each among two charts belonging to the same reference frame, 
and in each of which the metric has the form ([3]), then the two tetrad fields 
obtained thus are related together by a constant Lorentz transformation. It 
follows that, by choosing the diagonal tetrad in a chart, in which the met- 
ric has the form ([3]), we solve the non- uniqueness problem in any reference 
frame. In Sectional we review the procedure and its uniqueness features and 
we show the physical relevance of the metric (12]) • Appendix |X] explains our 
restriction to first-quantized Dirac theory. 



2 Conditions for equivalent operators 

For a wave equation on a spacetime manifold, it is the local expression of the 
wave equation which can be rewritten in the Schrodinger form: 

<9\I> 

(4) 

to get the Hamiltonian operator H that does not involve time derivatives [13J. 
Specifically, the Dirac wave function ip is a section of the spinor bundle E. 
The local expression of the Dirac operator involves choosing a local coordi- 
nate system on V: X i— >• (a; M ), and a local frame field (e a ) on E: X \-t (e a (X)). 
Thus, in Eq. (pfl), t = x°/c is the coordinate time and \I/ = (^ a ) a =i,...,4 is a 
function taking its values in C 4 , made with the components \I/ a of ip in the 
frame field (e a ) on E. We assume that there exists a global orthonormal 



3 The form ([3]) is well known |22) . It has been investigated previously in the context 
of the Dirac equation (23J El], but this was in the case of a stationary metric, g^,o = 
(thus here / t o = h t o = 0). For a stationary metric, independently of its form, Eq. (TIJ 
makes it natural to choose time-independent Dirac matrices. With this natural choice 
for a stationary metric, the non-uniqueness problem discussed here does not exist for the 
DFW equation 
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tetrad field on the tangent bundle TV. Then both the trivial bundle V x C 4 
and the complexified tangent bundle TcV are spinor bundles [BJ. Thus we 
have two different choices for E: i) we can take E = V x C 4 . This is the 
"quadruplet representation of the Dirac field" (QRD) j6j, for which the wave 
function is a complex four-scalar field. This choice is appropriate, in par- 
ticular, for DFW [BJ. (The DFW wave function transforms as a scalar on 
a change of the coordinate system [31 El ES]-) Or instead ii), we can take 
E = TcV. This is the "tensor representation of the Dirac fields" (TRD), for 
which the wave function is a complex vector field [6j [7J. 

The operator H depends in general on the coordinate system or chart, but 
(for a given 7 field, see Note [TJ), it remains the same for two charts which 
exchange by a purely spatial change [HI IT3] : 

x '° = x °, .,-'< piU'l) (j, k = 1,2,3). (5) 

The relation ([5J between two charts is an equivalence relation for charts which 
are all defined on a given open set U C V. We formally define a reference 
frame as an equivalence class F of charts for this relation, for a given open 
set U (the domain of F) [17]. This definition corresponds with the notion of 
a reference frame as being a fictitious fluid: the world lines 

x° variable, x J = constant for j = 1,2,3 (6) 

are the trajectories of the particles constituting the reference fluid [IB] . Our 
definition assumes, moreover, that the time coordinate: X h> x°, is a fixed 
function. This is made necessary by the fact that the Hamiltonian operator 
depends on the choice of the time coordinate. We emphasize that the de- 
pendence of the operator H on the reference frame is necessary and natural: 
it occurs for any wave equation [13]. That dependence has thus nothing to 
do with the dependence of H on the coefficient field 7^, Q which is the non- 

4 For the alternative versions of the generally-covariant Dirac equation [BJ, the r.h.s. is 
augmented with a term that involves the hermitizing matrix field A [U [2l 116] : 

Y l D^ = - ^A-^D^Af))^, (7) 

in order to ensure the current conservation. That term vanishes for the DFW equation 
[16] . The matrix A is determined by the data of 7^ only up to a positive scalar field X(X) 
[16) . Moreover, A is crucial in the definition of the scalar product [16j |9], also for DFW 
with a general set (7 t,Q ) in Eq. (|29|) . Thus, A is an additional coefficient field. Upon a 
local gauge transformation S, the field A changes to A = S^AS [21 Ej. 
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uniqueness problem that occurs specifically for the generally-covariant Dirac 
equation, and which is related to the fact that it is a gauge theory. 

As shown in detail in Ref. [6J, the connection D used to define the Dirac 
operator can be any connection on the spinor bundle E (for which there 
are two different choices, see above). The general expression of the Dirac 
Hamiltonian [5] involves the matrices V ^ of the connection D and/or the 
corresponding covariant derivatives acting on = <9 M + r M . On changing 

the coefficient fields 7 M [and A, Note H] by a local gauge transformation S 
according to Eq. there are two distinct possibilities for the matrices 

• Either [let us call this "Choice (i)"] they change according to the rule 

r^s^r^s + s-^s), (8) 

which is the necessary and sufficient condition in order that the Dirac 
equation be covariant under the gauge transformation fl2j), accompanied 
by the change 

* = S~ l y (9) 

of the wave function. This rule applies automatically if the transfor- 
mation d2J), © is "passive", i.e., if it results from a mere change of the 
frame field on the spinor bundle E, the connection being left unchanged 
[BJ. On the other hand, consider an "active" gauge transformation, i.e., 
one for which Eqs. (jSJ) and (jUJ) relate the expressions of, respectively, 
two different Dirac 7 M fields and two different Dirac ip fields, the frame 
fields on E and (for 7^) on its dual E° and on TV being left unchanged 
[BJ. In that case, the rule (|SJ) means that the connection D on E is 
changed for another connection D on E. The rule (|5D is verified for 
DFW by the very construction of DFW [8], whether the gauge trans- 
formation fl2}, (J9D is passive or active. 

• Or still ["Choice (ii)"], the matrices may be left invariant [TBI [9], 
meaning that the connection on the spinor bundle E is left unchanged, 
the frame field being also left unchanged: 

f „ = rv (10) 

Although this condition does not generally ensure the covariance of 
the Dirac equation under the gauge transformation (j2J), (JSj) [IB], it 
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does make sense: any two sets of Dirac matrices fields, say 7^ and 
7 M , both satisfying the anticommutation relation (TjQ), are equally valid 
coefficient fields for the Dirac equation and thus for the Hamiltonian 
derived therefrom, whether or not they lead to equivalent Dirac equa- 
tions. Choice (ii) is appropriate to investigate the alternative versions 
of the curved-spacetime Dirac equations, because each of them is based 
on using a specific connection on the relevant bundle E [HI CTJ. Thus, 
each of them has the property that, in contrast with DFW, it does not 
lead to a unique Dirac equation in a given spacetime (V, g) [16j [6] , so 
that Choice (i) does not make much sense for these alternative versions. 

In Ref. [9], the conditions on the gauge transformation S have been 
derived in order that the Hamiltonian H got after the application of 5" be 
physically equivalent to the starting one, i.e., in order that 

B. = S~ 1 YLS. (11) 

This has been done separately for, on one hand, DFW, giving the very simple 
condition 

d S = (DFW), (12) 

and for, on the other hand, two alternative versions of the Dirac equation 
[TJ, based on TRD. It is easy to check that the derivation made for the TRD 
case, and thus the condition derived for that case [9J, applies in fact also, as 
it is, to "Choice (ii)" above, independently of whether QRD or TRD is being 
considered. Thus, for Choice (ii), we have the following condition in order 
that Eq. ffTTj) be satisfied: 

B°(d S)S- 1 = [B^D^S- 1 ]* (Choice (ii)), (13) 

where B^ = A'~f fl , A being the hermitizing matrix field [H [21 [16], and where 
M a = |(M-M f ) is the antihermitian part of a matrix M, with M f = (M*) T 
denoting its Hermitian conjugate. Moreover, D^S is the 4x4 matrix made 
with the covariant derivatives D fJi S a b of the tensor field 

S = S\e a ®e\ (14) 

where (e a ) is the chosen frame field on E and (6 b ) is the dual frame field on 
the dual vector bundle E°, and is given by [6] 

d^s = (d^s\) = d,s + r^s - s r M . (15) 
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Note that, from (1141) . it follows that the complex matrix defining a local 
similarity transformation, S = {S a b ), behaves as a scalar under a coordinate 
change, the frame field (e a ) on E being left unchanged. 

In the case of time-dependent fields, the Dirac Hamiltonian H is in general 
not Hermitian [THIEZ]- Then, the relevant spectrum is that of the Hermitian 
part of H: 

E = H s ^i(H + Ht), (16) 

where the Hermitian conjugate operator H* is with respect to the unique 
scalar product identified in Ref. [16]. The name "energy operator" for the 
operator E is justified, see Ref. [9] and Eq. ( IT6]) below. As with Eqs. (jl~2l) 
and (fl~3l) for H, in Ref. [9], the conditions on the gauge transformation S have 
been derived in order that the energy operator E got after the application of 
S be physically equivalent to the starting one, i.e., in order that E = S 1-1 E S. 
These are: 

[B°(d S)S- 1 ] a = (DFW), (17) 

and 

[B^D^S- 1 - B°(d S)S- l ] a = (Choice (ii)). (18) 

Hence, for DFW, both the Hamiltonian and the energy operator stay 
invariant in the sense of Eq. (TTTT) after a gauge transformation ([2]) applied to 
the field 7 M , provided that the gauge transformation S is time- independent, 
Eqs. (jl~2l) and (TTT|) . However, the validity of the condition doS = depends 
on the coordinate system. The chain rule shows that, in order to ensure that 
it applies in any chart, the stronger condition 

d^S = = 0, ..,3) (19) 

has to be valid. Thus, for DFW, the wish to get equivalent Hamiltonian 
operators and equivalent energy operators before and after a gauge transfor- 
mation, independently of the chart, leads us naturally to impose the restric- 
tion that the gauge transformation should be a "global" one, i.e., S constant 
over the spacetime. On the other hand, for Choice (ii), the presence of the 
covariant derivatives D^S in the corresponding Eqs. ffTSl) and (fT8"|) makes 
the situation more complex. However, there is one alternative version of the 
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generally-covariant Dirac equation for which it is simple. This is the "QRD 
0" version [5], which is got by taking the trivial connection on the bundle 
E = V x C 4 : 

= in the canonical frame field {E a ), (20) 

where the canonical global frame field X i— > (e a (X)) = (E a ) on V x C 4 is the 
constant canonical basis (E a ) of C 4 , i.e., E a = (5^). Thus, for QRD-0 with 
the canonical frame field, we have D^S = d^S from (1TB} and (12"0~1) . Hence, in 
that case, it is obvious from Eqs. (fT3"j) and (TIB"]) that d^S = is a sufficient 
condition in order that the gauge transformation S lead to a Hamiltonian 
operator equivalent to the starting one, and the same is true for the energy 
operator. 

Therefore, we will try to define a natural class of coefficient fields 7 M [and 
A, NoteH], any two of which exchange by a constant gauge transformation S. 
If we succeed in that, then, in any possible reference frame, one and the same 
operator H [in the sense of Eq. (TlT]) ] will be defined from any two choices 
of coefficient fields 7^ and 7^ taken from that class, as well as one and the 
same operator E. Moreover, this will be true for both DFW and QRD-0. 

3 Tetrads come from a square root of the 
metric 

In a given chart x '■ X a tetrad field (u a ) is given by its matrix 

a = (a^Q,), such that u a = a^ a d^. The orthonormality condition for (u a ) 
writes then 

g(u a ,up) = g^a^ a a v p = 7] aP , i.e. a T Ga = r], (21) 

where rj = diag(l, —1, —1, —1) is the Minkowski metric and G = (g^) is the 
metric's matrix in the chart x- Most of the explicit constructions of a tetrad 
field in a given spacetime concern the case of a "diagonal metric", i.e., one 
whose matrix is diagonal in the chart considered: 

G = diag(d M ), (22) 

that chart being admissible [18], i.e., g 00 > and the 3x3 matrix (gj k ) (j, k = 
1,2,3) is negative definite, thus here d > and dj < (J — 1,2,3). The 
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construction then consists (e.g. Ref. [28]) in defining the tetrad from the 
following matrix: 

a=(a^)=diag(l/ v /^). (23) 
This defines the "diagonal tetrad" 

u a = 8^djyj\d^\. (24) 

That particular tetrad is orthonormal for the particular metric ([22]) . because 
it satisfies the general condition (I2TI) . 

Now consider a chart in which the metric's matrix G has a general form, 
though satisfying the admissibility condition for the chart, just recalled. In 
that general case: i) as for the standard Cholesky factorization, one may show 
with the Minkowski metric 7] (replacing the identity matrix of the standard 
Cholesky method) that there exists one and only one solution of the equation 

b T v b = G (G — G(X),X e V), (25) 

which is a lower triangular matrix b = C (X) with (strictly) positive diagonal 
entries, ii) Any other solution b of fl25|) has necessarily the form 

b = LC (26) 

with L = L(X) some Lorentz transformation, L(X) £ 0(1, 3). 

Note that any solution b of fl25l) may be termed a "square root of the 
metric" , since b ■ c = b Tr qc is bilinear. These solutions are in one-to-one 
correspondence with orthonormal tetrads (u a ) in a metric whose matrix is 
the general G, because Eq. f[2B]) for b is equivalent to the condition (I2T]) for 

a = b~\ (27) 

The Dirac 7 M field associated with a tetrad field is (see e.g. Refs. [3J, [SJ): 

-f(X) = a» a (X) t*", (28) 

5 These two precise points were stated to the author by Frank Reiner [29]. Point (ii) 
is easy to check: if b and b' verify l[25|) . then L = b'b^ 1 verifies L T r] L = r\. Point (i) is 
proved in Ref. [50] , 

6 Works by Ogievetskii & Polubarinov i3T] and Pitts [35] are relevant in this connection 
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where the constant matrices y^ a (a = 0,...,3) obey the anticommutation 
relation with the Minkowski metric: 

7 ii« 7 W» + 7 ^ 7 Ha = 2ri af) U, a, /3 E {0, 3}. (29) 

Thus, for a general form of the metric, we are able to define a field of Dirac 
matrices, by solving Eq. fl25|) and then by defining the matrix a of an or- 
thonormal tetrad field according to Eq. (|27p . Furthermore, we seem to have 
a preferred choice for b in Eq. ff27|) : namely, b = C(X), the unique Cholesky 
decomposition of the metric's matrix G(X). Since a diagonal matrix is a 
fortiori triangular, C = diag(v/|d M |) is the Cholesky decomposition of the 
metric (|22|) . Hence the "diagonal tetrad" fflty associated with a diagonal 
metric, which in our opinion is a natural choice for that case, is a particular 
case of this choice b = C. But, in the general case of a metric that is not 
diagonal in the chart considered, there does not seem to be any physical 
reason (as opposed to a merely computational one) to prefer the tetrad field 
associated with C(X) over any other possible tetrad field. Even apart from 
this, there is a serious difficulty, which we analyze in the next section. 

4 Non-uniqueness of the "Cholesky prescrip- 
tion" 

What is physically given is the reference frame (a three-dimensional congru- 
ence of time-like world lines), not the coordinate system, for which there is 
a vast functional space of different choices within a given reference frame. 
As we recalled in Sect. [21 the Hamiltonian and energy operators H and E, 
hence the energy spectrum, depend naturally (for any wave equation) on the 
reference frame as we define it. That is, H and E depend naturally on the 
equivalence class modulo the purely spatial coordinate changes (jSJ). Unfor- 
tunately, for the curved-spacetime Dirac equations, the operator H (as well 
as the operator E) is not unique even in a given chart, hence it is a fortiori 
non-unique in a given reference frame. What we are trying to do, is to find 
a particular prescription — a particular way of choosing the gamma field in 
a given chart — ensuring the physical uniqueness of each of these operators. 
Thus, if we have two admissible charts x an d x' that belong to the same 
reference frame, applying a such prescription successively in the chart x an d 
in the chart x' should provide operators H and E, then H' and E', with H' 
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(respectively E') being physically equivalent to H (respectively E). Setting 

P^^TTT, (30) 

we have from ([5]): 

P° = l, P° i = P^ = 0(j = l,2,3). (31) 

Let us calculate the local transformation L that transforms a first tetrad 
(u a ), having matrix a in the chart x, to a second one (n«), having matrix 
a' in the chart x' — so that u'p = L a /i u a . From u'p = a'^gd^, d' u = P^ v d^ 
and <9 M = b a ^u a (where b = a" 1 ), we find immediately [independently of the 
validity of Eqs. © and §3)]: 

L% = b\P^a% } or L = b P a' . (32) 

If both tetrads are orthonormal, L given by ( 132]) is automatically a Lorentz 
transformation. In a general spacetime, b and a' depend on the common 
time coordinate x° = x'° as do G and G', while P does not depend on it by 
virtue of It is a priori unlikely that the time-dependences of b and a' 
cancel one another in Eq. ( 1321) . so we expect that in general L depends on x°. 



Let us begin to check this when, more specifically, each of the two tetrads 
(u a ) and (u'a) is the (orthonormal) "Cholesky tetrad" in the chart x (respec- 
tively x')- Thus, a = C^ 1 and a' = C'~ x , where C (respectively C) is the 
unique lower triangular matrix with (strictly) positive diagonal entries such 
that C T r]C = G (respectively C' T r]C' = G'), G and G' being the metric's 
matrices in the charts x an d x'> respectively. Note first that a and a' are 
lower triangular matrices, as are C and C. That is, a^ a = if fi < a. It 
follows that a M M = l/C^ at fixed \i (without sum), so that a (and a') has 
(strictly) positive diagonal entries, as has C (and C). In Eq. (1321) . we have 
now C = a -1 in the place of b. Since C and a' are lower triangular matrices, 
( |32l) rewrites as 

a 3 

By using (I3TI) . this specializes, on one hand, to 

3 

L % = C °o E a " ^ L °o = C °o a '°o and L° p = (p = 1, 2, 3), (34) 

v=(3 
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thus also (L _1 )° p = (p = 1, 2, 3), from which, using L T = r}L~ x rf x , we find 
L p = 0. Entering this into L T i]L = 77 gives us {L° ) 2 = 1, hence, in view of 
L° = C° a'°Q > 0, L° = 1. Thus we have simply: 

L° = l and L° p = L p = (p= 1,2,3). (35) 

On the other hand, for the purely spatial components of L, f )33|) becomes, 
using again (I3"T|) : 

P 3 

iP ff = E C7P *E P * a,i « (P,ff=l,2,3). (36) 

fc=l j=q 

In particular, for g = 3, we get 

L P 3 = a%irC\P k 3 (p= 1,2,3). (37) 
fc=i 

Moreover, since C is lower triangular and verifies C T r]C = G, we have 

3 

E^i = (7,* = 1,2, 3), (38) 

p=j 

so that c T c = —g, with c = (C P j), is the standard Cholesky decomposition 
of the positive definite symmetric 3x3 matrix h = (—gjk)- The same applies 
with primes, of course. In particular: 



1 

a 



— = C'\ = = V -P\ P\ 9jk > 0. (39) 



(Of course, unless explicitly mentioned otherwise, there is summation over 
repeated indices, here the spatial coordinate indices j and k varying from 1 
to 3.) Hence from (l3j 



V^P fiP pk 

L\= ^fc=i _fc£j_ (p= 1,2,3). (40) 

-P\ P\ g jk 

Clearly, this depends in general on the time coordinate as do the C p k 's and 
the g jk 's. 
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Even more specifically, let us assume that the metric has the diagonal 
form (1221) in the chart x (but not necessarily in the chart x')- Thus C = 
diag(yjd^l), so Eq. ( 140]) becomes 



L p 3 = J- ^ y^> (no sum on p= 1,2,3). (41) 
-{P\fdj 

Since P does not depend on x°, we find easily from this that 

° (L p 3 ) cc PMP^) 2 ^ (%) (no sum on p = 1, 2, 3), (42) 



<9x° v ^ ^ 6/ dx° \d p/ 

with a non-zero proportionality factor. In general, the eigenvalues dj, and 
their ratios as well, depend on x°, hence -£p (L p 3 ) will be non-zero. Two nat- 
ural exceptions are: i) dj = ojj(x)(^(x ) [x = (x J )] — which seems however 
difficult to get — and ii) dj = d°j h(x°, x) with d°- constant (d° < with h > 0). 
In the latter case, we get d'j = —h (j = 1, 2, 3) after the coordinate change 

x'- 7 = yj—dj. That is, Case (ii) corresponds precisely with the metric ([3]). 

Thus, in a non- stationary spacetime, the local Lorentz transformation L 
which transforms the Cholesky tetrad (u a ) got in a chart x to the Cholesky 
tetrad (u' a ) got in a chart x' that belongs to the same reference frame, in 
general depends on the time coordinate. [L is given by Eqs. ( 1315]) and (|36p . 
and this dependence is particularly obvious on Eq. (|4T)|) and its specialization 
for the case that G is diagonal, Eq. ( 142]) .] The fields of Dirac matrices 7^ 
and 7^ associated respectively with the two tetrad fields (u a ) and (u' a ), Eq. 
( 128|) . are related together by Eq. ([2]), where S is a local gauge transformation 
obtained by "lifting" the local Lorentz transformation L. That is, S is such 
that 

(VX G V) A(S(X)) = L(X), (43) 

where A : Spin(l, 3) — > SO(l, 3) is the two-to-one covering map of the special 
Lorentz group by the spin group. In other words, S is such that 

S(X) = ±S(L(X)), (44) 

where L 1— > ±S(L) is the spinor representation (which is defined only up to 
a sign, contrary to the covering map of which it is the "inverse") [16] . (We 
assume that a smooth lifting is possible, i.e., that 5* depends smoothly on 
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X G V; this is always possible in a simply-connected spacetime [5].) Since, 
as we have shown, the local Lorentz transformation L in general depends on 
the time coordinate, so also does the local gauge transformation flU]) . It fol- 
lows that the Dirac Hamiltonian operators H and H associated, in the chart 
X, with the fields of Dirac matrices 7^ and 7 M , are in general not physically 
equivalent — at least for DFW, see Eq. (I12p . However, both H and H are 
got by applying the "Cholesky prescription" a = C~ x in two different charts 
belonging to the same reference frame. Therefore, that prescription is not 
unique. 



5 Case of a space-isotropic diagonal metric 

As we have just seen, the "Cholesky prescription" a = C^ 1 does not in 
general provide a unique Dirac Hamiltonian operator H (nor a unique energy 
operator E) in a given reference frame, even if the metric is diagonal — in 
which case that prescription coincides with the "diagonal tetrad" prescription 
(}24j) . which is thus non- unique in general. Therefore, we are led to consider 
a special form of a diagonal metric. For the reason that appeared after Eq. 
(|4T|) . we are in fact led to consider the metric ([3]). And indeed, for this metric, 
the "diagonal tetrad" prescription is unique: 

Theorem 1. Let the metric have the space-isotropic diagonal form |3j] in 
some admissible chart x '■ X 1— > Let be an admissible 

chart that defines the same reference frame R, i.e., x and x' are defined on 
the same domain U and the change from (x^) to is purely spatial, Eq. 

(i) In order that the metric have the same form (TJjj in the new chart, it 
is necessary and sufficient that the transformation of the spatial coordinates 
be a constant rotation, followed (or preceded) by a constant homothecy. 

(ii) If one applies the "diagonal tetrad" prescription (Efy in each of the 
two charts, the two tetrads obtained thus are related together by a constant 
Lorentz transformation L, hence give rise, in any reference frame F, to equiv- 
alent Hamiltonian operators as well to equivalent energy operators — for the 
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DFW and for the QRD-0 versions of the Dirac equation. 



Proof. Under the assumptions of the theorem, the matrix P, with compo- 

dx^ 

Q x tv 1 



nents P^ v = f 2 ^, has the form (13T1) . thus it has the block structure 



p= {oq)' < 45 > 

where Q is the 3x3 matrix with components Q 3 k = J^. In the chart 
the metric has the space-isotropic diagonal form ([3]), thus in particular the 
metric's matrix G has the same block structure as has P. Moreover, for G, 
the 3x3 matrix which corresponds to Q in Eq. ( l4"5j) is just — /1I3. Hence, 
the metric's matrix G' in the new chart still has that block structure and 
writes explicitly: 

G' = P T GP=(^ _ h ° QTQ ). (46) 

It thus has the form ([3]), iff there is some function h' = h'((x'^)) > [the 
new coefficient for the spatial part of the metric (J3J)], such that 

hQ T Q = h'l 3 , (47) 



or 

C^F-F^U, F = Q" 1 , F> t= °*. (48, 

Recall that the coordinate change is assumed purely spatial, so F and C do 
not depend on the time coordinate x°. Hence, for instance, even though h 
depends in general on x° as well as on the spatial coordinates, the ratio h/h' 
cannot depend on x° if Eq. (I48p is verified. 

From now on, the proof of point (i) is based on a mathematical re- 
sult that pertains to continuum mechanics, namely it uses the theory of 
the compatibility of a field of large deformations. Consider a continuous 
medium whose "initial configuration" is described by some open domain Q 
of R 3 , and let C be a candidate field for the "right Cauchy- Green strain 
tensor" of that continuous medium — thus C is a symmetric, positive defi- 
nite, spatial tensor field of the (0 2) type, defined over the domain Q. The 
compatibility problem is the question whether there exists a transformation 
(f : Q — > R 3 , x = (x 3 ) h-> x' = (x' J ), which admit this tensor field C as 
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its Cauchy-Green strain tensor, i.e., which be such that the gradient ten- 
sor F of the transformation, with F\ = verify indeed F T F — C. A 
tractable answer to this question has been given by Vallee and collaborators 
I3U [35]. This answer is as follows: 



a) Let U = U(x) be the unique positive square root of C, i.e., U is 
symmetric, positive definite, and satisfies U 2 = C. If there exists a trans- 
formation ip whose gradient F verify F T F = C, set R = FU~ l . We get 
R T R = 1 3 , thus R is an orthogonal transformation: R = R{x) G 0(3). 
That is, 

F = RU (49) 

is the (unique) "left polar decomposition" of F. The spatial variation 5R of 
R in an arbitrary direction <5x provides an antisymmetric tensor field: 

R T 5R = -(R T 5R) T } (50) 

and therefore has the form 

R T 5R = j(A(x)6x), (51) 

where v h-> j(v) is the linear correspondence which associates, with a vector 
v of R 3 , the antisymmetric tensor T such that Vx G R 3 Tx = v A x, with 
A denoting the vector product. Thus, x 1— > A(x) is a field of linear transfor- 
mations of R 3 . A necessary condition for the existence of the transformation 
ip : x 1 — y x' is that A should be given by 



.4 



detU 



U(RotU)U- tl(UR ° tU >U 



(52) 



where Rot U is the field of linear transformations defined thus (at any given 
point x G Q): for any constant vector v G R 3 , 

(Rot U)v = rot (Uv) = V A (Uv). (53) 

b) Then, A being given by (|52|) . the necessary and sufficient condition for 
the local existence of the transformation ip is that (at any point x G fl): 

(Rot A T ) = -Adj A, (54) 



7 In continuum mechanics, R must be specifically a rotation, but this does not play 
any role in the derivation [33l [3H [35] . 
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where the "adjunct" linear transformation Adj A of the linear transformation 
A (or rather, of the linear transformation A(x) for a given point x) is the 
linear transformation B such that 

Vu, v,weR 3 (u,v,5w) = (iu,Av,w), (55) 

(u, v, w) = (u A v).w being the scalar triple product of three vectors of R 3 . 

Let us apply this theory to the case that C is given by Eq. ([15)1 . We have 
thus 

U = al 3 , a = (h/h') l l 2 . (56) 
We get from this and the definition f[53l) : 

(Rot U)w = V A (av) = (Va) A v, (57) 

(U Rot U)v = a (Va) A v, (58) 
thus Rot U and [7 Rot U are antisymmetric, and we have 

tr (17 Rot 17) = 0, det(£/) = a 3 . (59) 

Hence, from Eq. ([52)) . A is simply 

A = U- 1 Rot £7, Av = a' 1 (Va) A v, (60) 

and is thus antisymmetric. It follows that, for any constant vector w 6 R 3 , 

(Rot A T )w = -V A [a' 1 (Va) A w] 

= -(V.w) a" 1 Va + [V.(a _1 Va)] w 

= [-a' 2 (Va) 2 + a^Aa] w. (61) 

Let us compute £> = Adj A. Setting g = Va, we get from f[6T))) and the 
definition ([55)) : 

Vu, v,wGR 3 (u,v,5w)a 2 = (g A u, g A v, w) 

= [(g Au) A (g A v)] .w 

= {[(g Au).v]g- [(g Au).g]v}.w 

= [(g Au).v] (g.w) 

= (u,v,g)(g.w), (62) 
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that is, 

Bw = a" 2 (g.w)gE«" 2 [(Va).w]Va. (63) 
From (IB"Tj) and (I6"3"j) . the compatibility condition (I5"4j) is satisfied, iff 

Vw G R 3 , [-^ 2 (Va) 2 + a"'Aa] w = -a" 2 [(Va).w] Va. (64) 

Clearly, this happens if, and only if 

Va = 0. (65) 

If Va = 0, we have A = from (IBTfl) . hence from (15TI) : 

R = Constant. (66) 

It follows then from (j4*9l) and (156]) that 

F = = ai?. (67) 

<7X 

Note that, if it happens that deti? = —1, we may rewrite this as ^ = 
(— a)(— R). From ( 1651) and the remark after Eq. ( 1481) . it follows that, if Eq. 
(148]) is verified, then a = (h/h') 1 ^ 2 does not depend on x° either. Then, Eq. 
( I67j) means that indeed the transformation : x i— > x' is a constant rotation, 
followed (or preceded) by a constant homothecy. Conversely, if that is the 
case, we may write Eq. ( IBTl) . whence it is clear that Eq. ( 1481) is verified. 
Point (i) of lTheorem H is proved. 



Point (ii) now follows easily. First note that, applying a mere homothecy 
to the spatial part of a coordinate system in which the metric is diagonal, 
one leaves the diagonal tetrad (124"]) unchanged: 

K = d \J = d J yK\ = u » ( no sum ) • ( 68 ) 

Hence, we may assume that there is no homothecy, i.e., a = 1 and R G SO(3) 
in Eq. ( 1671) . Thus F = R in Eq. ( piBl) . hence h' = h [so the metric ([3]) is left 
invariant by the coordinate transformation], and dx'i /dx k = R 1 k . Hence by 
(El: 

d -R J &■ 

Uk = " = — y=^- = R 3 \u'a (no sum on k = 1, 2, 3). (69) 
V-dk Vh 
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Together with (1531) . this means that one goes from the diagonal tetrad as- 
sociated with the first coordinate system to the diagonal tetrad associated 
with the second coordinate system by the constant Lorentz transformation 

fl 0\ 

The local gauge transformation (|44p obtained by "lifting" that constant L is 
of course constant: d^S = 0. (This is true in any chart, whether it belongs to 
the frame R or not, since, as we noted after Eq. ( IT5|) . S behaves as a scalar 
under a coordinate change.) As we concluded at the end of Section [2j this 
is a sufficient condition in order that the gauge transformation S lead to a 
Hamiltonian operator equivalent to the starting one, and also to an energy 
operator equivalent to the starting one, for both DFW and QRD-0. This 
completes the proof of lTheorem II 



6 Discussion 

i) The proposed prescription and its uniqueness. Let us assume that the 
metric has the space-isotropic diagonal form ([3]) in some coordinate system 
(chart) x '■ X i— > (x M ). [One may assume that x is a global chart, for the 
simplicity of discussion.] That chart defines a preferred reference frame, say 
R: this is the class of the charts that exchange with x by a purely spatial 
coordinate change ([5]). We define the field of Dirac matrices in the chart x 
from the diagonal tetrad field (u a ) defined in that same chart x- That is 
[Eqs. d23]) and (|28|)]: 



I ( no sum), d o = f, di = d 2 = d 3 = -h, (71) 

where (7^ a ) a =o,...,3 is a set of constant "flat" matrices, i.e., obeying the an- 
ticommutation relation ( |29l) with the Minkowski metric. In an arbitrary 
reference frame F, the field of Dirac matrices, say C 1 , is defined by trans- 
forming it as a four- vector [3j [6] from the chart x ^ R to an arbitrary chart 
cf) \ X i — y (y 11 ) belonging to the frame F: 

r=§£r- (72) 
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The two fields 7 M and C M are the expressions in the charts x an d 0, respec- 
tively, of a unique "7 field" (Noted]): the one associated with the tetrad field 
(u a ) {Ref. [6], Appendix}. [The fields (u a ) and 7 are global if the chart 
X is global.] Thus, for each choice of a chart x that both belongs to the 
reference frame R and is such that the metric has the form ([3]), we have one 
preferred tetrad field (u a ) and one preferred 7 field. Point (ii) of lTheorem II 
proves that the choice of a such chart x h as no influence on the Hamiltonian 
and energy operators. By the way, note that the invariance of the similarity 
matrix S under a coordinate change, used at the end of the proof of that 
theorem, is easily rechecked from (j2J) and ( 172]) . 

Of course, besides changing the chart in the particular reference frame 
R, one may also change the chart in the general reference frame F, say from 
4> : X 1— > (?/') to 0' '. X 1 — y But this does not change the 7 field. In 

other words, the two fields of Dirac matrices and ( >fl got by Eq. (1T21 are 
related together by the 4-vector transformation: C,'^ 1 = fj^rC \ with a trans- 
formation 1 — y (y ffl ) that is purely spatial — so that, as recalled in Section 
El the Hamiltonian operator, as well as the energy operator, remain the same 

The choice of the set of constant "flat" matrices, i.e. the set (7^ a ) a =o,...,3 
in Eq. ( 128]) . has no influence either. Indeed, two different choices, say (7^) 
and (7"°), give by ( 128]) two fields of Dirac matrices which exchange by the 
same constant similarity that transforms (7^) to (7^). 

All of these results are valid for both the standard version of the curved- 
spacetime Dirac equation ("DFW") and the alternative version with zero 
connection matrices ("QRD-0"). These two versions are in general not equiv- 
alent: although the DFW equation is equivalent to a particular case of the 
QRD-0 equation [26], this equivalence takes place after changing the field of 
Dirac matrices in a way that has no reason to be compatible with our pre- 
scription. Thus, our prescription leaves us with two non-equivalent choices for 
the Dirac equation in a given spacetime (V, g), leading to two non-equivalent 
choices for the Dirac Hamiltonian and for the Dirac energy operator in a given 
reference frame. Experiment might decide between the two choices. 

8 Two points may further be noted in this connection, i) Not only a unique DFW 
equation (as is usually the case) but also a unique QRD-0 equation is got by the proposed 
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ii) Physical relevance of the metric (TJ)). The spacetime metric cannot gener- 
ically be put in the form Q in a coordinate chart for general relativity. 
Indeed, even the standard post-Newtonian (PN) metric, which is valid for 
the usual, weak gravitational fields, has non-diagonal terms: specifically it 
has terms g$j (j = 1,2,3). However, apart from these terms it has just the 
form 02]), moreover the g$j terms are 0(c~ 3 ) and are very small as compared 
with the leading corrections to the Minkowski metric, which are 0(c~ 2 ) and 
depend on the Newtonian potential [36] [37] . The energy levels of a Dirac 
particle, as obtained with the 0(c~ 2 ) corrections to the Minkowski metric, 
differ already very little from the energy levels got with the non-relativistic 
Schrodinger equation in the Newtonian potential [23] EH]- In addition, the 
metric ([3]) , assumed valid in the reference frame R, transforms to a more gen- 
eral form with non-zero goj terms, after a boost and/or a rotational motion. 
If these involve linear velocities of post-Newtonian order: v = O^ 1 ), then 
one gets g^j = 0(c~ 3 ) as in the standard PN metric. Thus, the form ([3]) is 
general enough for the prospective purpose of testing the generally-covariant 
Dirac equations in a realistic spacetime metric. ^ 

iii) More general metrics and other prescriptions ? The metric ([3]) is indeed 
valid in a preferred reference frame in an alternative, scalar theory of grav- 
ity, and this in the most general case j3H]- Within general relativity (GR), 
on the other hand, one may consider it frustrating that a unique prescrip- 

prescription. Indeed, two different choices of the chart x € R lead to 7 fields exchanging 
by a constant similarity transformation S: d^S = 0. Now, since the QRD-0 equation fixes 
the connection matrices (to zero), the condition of covariance of the Dirac equation 
under a similarity is derived in the same way as in Ref. |16j . and is thus the same {Eq. 
(69) in Ref. [16]}, i.e.: S^B^iD^S) = (D^B^S, where B» = . Thus, with r„ = 0, 
this is verified by S since d^S = 0. ii) Any form of the 4-scalar Dirac equation ("QRD") is 
equivalent to a 4- vector Dirac equation ( "TRD" ) , essentially by switching from the global 
frame field (E a ) on V x C 4 to an (arbitrary) global frame field (e a ) on TcV ([6i, Theorem 
1). This TRD equation, and the Hamiltonian and energy operators associated with it 
in a given reference frame, have identical local expressions with those of, respectively, 
the starting QRD equation and the corresponding H and E operators. Therefore, the 
two non-equivalent choices for the Dirac equation in a given spacetime, got by using our 
prescription, can be regarded as pertaining to TRD as well as to QRD. In short, we have 
two choices only, regardless of the four-scalar or four- vector representation. 

9 See Ref. [30] for a proof that our assumption: (A) "There is a local chart x : V D 
U — > R 4 , such that the local expression of the metric g is ((3]) in that chart," is invariant 
under (isometric) diffeomorphisms. 
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tion has been found only for that metric ([3]) — even though, as shown above, 
this metric is general enough for the foreseeable experimental tests. How- 
ever, recall that we were naturally led to this metric in Section H] in order to 
get unique Hamiltonian and energy operators from choosing the "diagonal 
tetrad". Indeed, we do not see any other general solution to get (L p 3 ) = 
in Eq. (|42p . Moreover, to have this is just a necessary condition in order 
that the Cholesky prescription be unique for some diagonal metric. For the 
well-known exact solutions of GR, either the metric can be put into the form 
as is the case for the Schwarzschild spacetime (choosing isotropic coor- 
dinates) and the Friedmann-Lemaitre-Robertson- Walker spacetimes, or it is 
unlikely that the Cholesky prescription is unique, as for e.g. the Kerr space- 
time. One may also examine other prescriptions. There is indeed one known 
other prescription which provides a unique tetrad field in a given coordinate 
chart. This is the "spatially symmetric time gauge" [HUE]) which imposes 
to the tetrad's matrix a two conditions: (i) a° p = (p = 1,2,3) ("time 
gauge" [32], leading to what is called a Schwinger tetrad in Refs. [HI IS]; see 
also Ref. [13]); and (ii) the matrix (a J p ) (j,p = 1, 2, 3) is got from the (sym- 
metric) positive square root, say a = Vh, of the positive definite symmetric 
matrix h = {—Qjk) (j,k = 1,2,3). [Note that, for a diagonal metric, thus 
in particular for the metric ([3]), both the Cholesky tetrad and the "spatially 
symmetric time gauge" tetrad coincide with the diagonal tetrad (124"1).] How- 
ever, the tetrads got by applying that prescription in two admissible charts 
X and x' that belong to the same reference frame F are still related by Eq. 
d32]). From fl3j]) we get ti = Q T hQ, with Q = (P j k ), thus ti = R T UhU R } 
where Q = UR is the right polar decomposition. Since the eigenvectors of 
U and h do not coincide for a general matrix Q, there is no simple relation 
between y/h and Vti, and thus there is no simple relation between a and a'. 
It is hence very unlikely that the time-dependences of b = a^ 1 and a' may 
cancel one another in Eq. (132]) for a general-enough metric differing from fl3]) . 
So that L should depend on x°, which would mean that Kibble's prescription 
[4TU E] is not unique and can hardly be hoped to be made unique except for 
the metric ([3]), just like the Cholesky prescription. 

Acknowledgement. I am grateful to Frank Reifler for allowing me to 
use a private communication and for remarks about the present paper, to J. 
Brian Pitts for pointing out Ref. [31] . and to both for exchanges about the 
non-uniqueness problem. 



24 



A Appendix: First vs. second quantization 



The non-uniqueness problem at hand has been proved to be there for the 
Hamiltonian and energy operators associated with the covariant Dirac equa- 
tion, thus for first- quantized Dirac theory [9]. So the present proposal for a 
solution remains within first-quantized theory. There is a well-known argu- 
ment against first-quantized Dirac theory: due to the negative-energy states, 
the positive-energy levels — e.g., of an atomic electron — would be instable 
since the system would make transitions to negative energy levels. However, 
it is also well known that, in low-energy situations, to which the use of the 
first-quantized Dirac theory is limited, the transition between positive and 
negative energy states can safely be ignored and consideration can safely be 
restricted to positive-energy states. For instance, the classical discussion of 
the hydrogen-type atoms uses the Dirac positive-energy states of the electron 
in the nucleus's electromagnetic field. At low energies, one thus may add a 
selection rule that only the positive-energy states are allowed. The quantum 
mechanics of the Dirac equation makes a scientific literature of tens of thou- 
sands of papers and books, including some of the pillars of modern physics. 
The energy levels of specifically the DFW equation were discussed by well- 
known physicists such as Brill & Wheeler [3], Greiner et al. [44J, Hehl & Ni 
[12], and Ryder [TO], to name just a few. Indeed, the three effects on quan- 
tum particles in the gravitational field which have been observed so far — the 
COW effect [15], the Sagnac effect [16], and the quantization of the energy 
levels [UJ — are low-energy effects and are discussed in first-quantized Dirac 
theory. This is also true for the Mashhoon effect [TH [12] , although testing 
the latter is still an experimental challenge. Obviously, the non-uniqueness 
of the first-quantized Dirac Hamiltonian operator and the energy spectrum 
is a serious problem for the discussion of these effects [9]. Moreover, note 
that the non-uniqueness problem is there in any admissible coordinate sys- 
tem in essentially any spacetime [9] . Thus, it is there already in a Cartesian 
coordinate system in a Minkowski spacetime [13] — as soon as one uses the 
DFW equation with its gauge freedom, instead of using the original Dirac 
equation. It is easy to check that this holds true in the presence of an electro- 
magnetic field [30] . Therefore, the classical discussion of the hydrogen-type 
atoms, which is based on the quantum-mechanical Hamiltonian and energy 
spectrum, cannot be done if one uses the DFW equation with its gauge free- 
dom. 
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One may ask whether, instead of restricting the admissible gauge trans- 
formations as was done here, this non-uniqueness could be solved by going 
to a second-quantized theory. However, a central object in quantum field 
theory (QFT) is the quantum energy-momentum tensor. This should be de- 
rived from the "classical" energy-momentum tensor, i.e. in fact, that of the 
first-quantized Dirac theory. Even in that "classical" theory, the definition 
of the energy-momentum tensor is not fully clear. The canonical definition 
of that tensor: 

"• = 5(W*» + ^e(W-* £ (73) 

is univoque at least when the Lagrangian L is fixed, and is often used in the 
literature (e.g. [27JET]). The following Lagrangian [9] extends the standard 
DFW Lagrangian [3j [27] to the case of a general hermitizing matrix field A: 



- {D^) ] A^m + 2imtfL4tfl , (74) 

and the associated Euler-Lagrange equation is indeed the general Dirac equa- 
tion ©. From dZ3J and (JZ2J, we get 

&Af(d v V) - (d v V)* A'fv] - 5%L. (75) 
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Since and A transform as scalars and 7 M as a vector under a change of the 
chart, (175p defines a generally-covariant (1 1) tensor. One may replace t^ v 
by its symmetric part, but this is indifferent in what follows. The problem 
is that the "field energy" E of the "classical" Dirac field obeying Eq. (J7J), is 
equal to the expected value of the energy operator E [H] : 

E = J t° y^ d 3 x = (* | Etf) ee J *U 7 °(E*) y/^g d 3 x, (76) 

as was shown by Leclerc [27] in a less general setting. [Using the fact that L 
vanishes "on shell", this is easy to check from ( 1751) and the definition ( I16j) .] 
Thus, the non-uniqueness of the energy operator E translates to that of the 
"classical" energy-momentum tensor, indeed even to the space integral of the 
latter. 

On the other hand, there is also Hilbert's definition of the energy-momentum 
tensor, as the derivative of the Lagrangian density £ = y/^g L with respect 
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to variations of the metric. This leads (up to the symmetrization) to just 
the same Eq. (1751) . though with the partial derivatives replaced by covariant 
derivatives [3l|52]. It is easy to check that, in contrast to the canonical tensor 
[Eq. fl75|) as it is], that other tensor is gauge- invariant, thus unique. How- 
ever, Eq. fl76|) then does not hold true (except for the non-generic case that 
the two tensors concide). I.e., the gauge- invariant field energy thus obtained, 
say E' ^ E, is not equal to the mean value | E\l/) of the energy operator 
( fl6|) . and therefore is not relevant. (If one would define a different energy 
operator E' 7^ E from its mean values E', this would have no relation to 
the Hamiltonian H, moreover it would be unfeasible to calculate its energy 
spectrum.) Since the same gauge freedom should apply to the first and sec- 
ond quantized Dirac theories, and since the predictions of the first-quantized 
theory should be deducible from the second-quantized ones in the low energy 
limit, this means that the relevant quantum energy- momentum tensor should 
be got from "quantizing" the canonical tensor — and thus should share this 
non-uniqueness with its classical counterpart. 

Moreover, in the works on QFT of the Dirac field in a curved spacetime 
(e.g. [HI HU |50]), we do not see any outline of the computation of one of the 
already-measured or prospective effects mentioned above [TTJ [121 HSJ HHJ HZ], 
which would allow one to hope that this effect might be unambiguously 
predictable in that framework. Finally, as noted above, the non-uniqueness 
problem found for the covariant first-quantized Dirac theory affects also the 
case with an electromagnetic field in a flat spacetime, even in an inertial 
frame. We therefore feel that one should first build an unambiguous first- 
quantized theory, as we attempted to do here. 
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